AD-AOHl  446 


UNCLASSIFIED 


SCIENCE  applications  TnC  HERKEI.EY  CALIF  F/6  9/1 

calculation  of  The  EOUIVA.FNT  capacitance  of  a rib  NFAR  a SINOL— ETC(U) 
JUN  77  s COEN*  T K LIU*  F M TESCH*:  F29601-76-C-012S 


AFVL-TR-77-60 


1fil4l44e 


AU  j ADA041446 

DDCiiLLM 


AFWL-TR-77-60 


AFWL-TR- 

77-60 


I 


CALCULATION  OF  THE  EQUIVALENT 
CAPACITANCE  OF  A RIB  NEAR  A 
SINGLE-WIRE  TRANSMISSION  LINE 

Science  Applications,  Inc. 

Berkeley,  CA  94701 

June  1977 


Final  Report 

Approved  for  public  release;  distribution  unlimited. 


L/ 


17 


D D C 

JUL  11  '9H 

JllcisignrE 


n' 


AIR  FORCE  WEAPONS  LABORATORY 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base,  NM  87117 


This  final  report  was  prepared  by  Science  Applications,  Inc.,  Berkeley, 
California,  unde'"  Contract  F29601 -76-C-0125 , Job  Order  12099513,  with  the 
Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  ’lew  Mexico.  Captain 
Harrison  (ELF)  was  the  Laboratory  Project  Officer-in-Charge. 

When  LiS  Government  dravnngs,  specifications,  or  other  data  are  used  for 
any  purpose  other  than  a definitely  related  Government  procurement  operation, _ 
the  Gcvernm.ent  thereby  incurs  no  responsibi  lity  nor  any  obligation  whatsceve*'. 
and  the  fact  that  the  Government  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawincs,  specifications,  or  other  data  is  not  to  be  regarded 
by  implicat''on  or  other.'mse  as  in  any  manner  licensing  the  holder  or  any  other 
person  :r  c:rp>'cation  or  conveying  any  r^gnts  or  permission  to  manufacture, 
use,  or  sell  any  patenteo  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  reviev/ed  by  the  Information  Office  (01)  and  is 
releasable  to  the  'iational  Technical  Information  Service  (,NTIS).  At  NTIS,  it 
will  be  available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


MICHAEL  G.  HARRISON 
Caotain,  USAF 

Project  Officer  FOR  THE  COMMANDER 


AARON  3.  L0GGI’<S;/Lt  Colonel,  LSAF 
Chief,  Phenomenolony  ^ Technology  3r 


/ 


AMES  L. 


olonal . 


GPIGGS,  JR. 
LISAF 


Chief,  Electronics  Division 


DO  NOT  RETURN  THIS  COPY.  RETAIN  OR  DESTROY. 


1 


UNCLASSIFIED 


SEC'jwiT’y  Classification  of  tmis  page  '’Whmn  o«<« 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETINO  FORM 

1 REPO^NUMBER  '2  r.nVT  Ar.r.FSSION  NO. 

AFWLfTR-77-60  J ; 

3 ftgCiPiEN'^’S  catalog  number 
/ 

A ^i^wE  mnd  Subtitle) 

5 ”YPE  OF  REPO®'^  4 PERIOD  COVERED 

Final  Sep^t  , 

CALCULATION  OF  THE  ,EQUI VALENT^CAPACITANCE  OF  A 
RIB  ^EAR  A JINGLE-i^IRE  TRANSMISSION  . — - 

k'  PERFORMtSG  ORG  Rtf^OR*'  NUMBER 

Shinion/Coen^ 

Tom  K./Liu 
Frederick  M.^Tesche 

9 COnTRaC^  DR  grant  suMBER^ft 

• F29601-76-C-P125  / 

? PERFORMING  OpG’aniZ  ATION  \amE  AnO  ADDRESS 

Science  Applications,  Inc. 

POB  277  ^ 

Berkeley,  CA  94701 

\Z  program  ElcMEn  ' prDjEC*  "a^^ 

AREA  4 WORK  UN  ''  numbers 

64747F  — 

^T2^513  / -■  - ' 

II  COnTPOllinG  OFFICE  name  anO  ADDRESS 

Air  Force  Weapons  Laboratory  (ELP)  ,/ 

Kirtland  AFB,  NM  87117  ' 

12  REPORT  date 

, Jun^^i977 

’3  NUMBER  OF  RAGES 

32 

•A  mONI’CRiNG  agency  name  S AOCRESSFI/  dillerent  'rom  C mtrolhng  Office) 

1 

15  security  Class.  »*>/  report) 

Unclassified 

'5«  declassification  downgrading 
SCnEDuLE 

*6  Distribution  statement  ©/  tht$  Report^ 

Approved  for  public  release;  distribution  unlimited. 

’T  Distribution  statement  at  the  eteirsct  entered  in  Block  ZO.  1/  dtl/ereni  from  Report) 

1»  SoPPLEMENTARY  notes 

19  i<Cy  WORDS  Continue  on  reveeee  side  i/  neceeserv  end  /d*nrf/y  by  block  number) 

Transmission  Lines 
Single-Wire  Transmission  Lines 
Capacitance 

Electromagnetic  Pulse  (EMP)  Coupling  to  Transmission  Lines 

i 

abstract  'Continue  »n  reeerse  side  It  necessmry  end  identity  by  block  number) 

The  effect  of  a rib  lying  on  a ground  plane  on  a single  wire  transmission 
line  can  be  represented  by  a lumped  capacitance.  In  this  paper,  the  capaci- 
tance is  evaluated  by  solving  a set  of  coupled  integral  equations. 

Parametric  studies  are  carried  out  and  the  results  are  presented  in  graphi- 
cal form. 

DD  » JAN  "1  1473  COITION  OF  1 NOV  S5  IS  OBSOLETE 


UNCLASSIFIED 


security  classification  '^Mis  B»JC  »T.en  D.IA  Em.rmi  ^ 

//  \ -i  f i I ,v 


V)1  < I j 


li 

Y 


WT 


4 

I 


TABLE  OF  CONTEN'T’S 


Section 

I 

II 

III 

IV 

I. 

V 

VI 

REFERE'ICES 


INTRODUCTION 


FORMULATION  OF  COUPLED  INTEGR.AL  EQUATIONS 


APPROXIMATE  SOLUTION  OF  THE  COUPLED 
INTEGRAL  EQUATIONS 


HERMITE  POLYNOMIALS 


NUMERICAL  SOLUTIONS  OF  THE  COUPLED 
INTEGRAL  EQUATIONS 


NUMERICAL  RESULTS 


1 


Page 

3 


10 

19 

21 
2 5 
32 


LIST  OF  FIGURES 


Figure  Page 

1 Geometry  of  the  Problem  5 

2 Equivalent  Circuit  of  the  Septum  Discontinuity  6 

3 Distribution  of  normalized  excess  charge 

for  h/b  = 0.1  . Mote  that  ix,  is  an  even 
function  of  x . 27 

4 Distribution  of  normalized  excess  charge  for 

h/b  =0.4.  " 28 

5 Distribution  of  normalized  excess  c.harae  for 

h/b  = 0.7  . 29 

6 Capacitance  due  tc  septum  as  function  of 
septum  height.  C^  is  capacitance  per  unit 

length  of  unperturbed  wire.  30 

7 Capacitance  due  to  septum  as  function  of 

septum  heiaht.  31 


2 


SECTION  I 


INTRODUCTION 

As  described  in  a number  of  recent  reports  (ref.  1,2,3) 
the  calculation  of  EMP  eneray  propagation  within  an 
electrically  complex  system,  such  as  an  aircraft,  often 
employs  simple  transmission  line  theory.  Such  propagation 
models  usually  consist  of  one  or  more  uniform  trans- 
mission lines  with  discrete  loads  and  distributed  sources 
arising  from  the  incident  electromagnetic  fields.  In 
order  to  more  accurately  account  for  the  nonuniform 
surroundings  of  actual  transmission  lines  on  aircraft, 
a number  of  "canonical  problems"  have  been  suggested  in 
ref.  (4).  By  solving  such  problems,  it  is  possible  to 
obtain  estimates  of  the  effects  of  local  perturbations  of 
the  transmission  line  fields,  and  later  relate  these  to 


1.  Tesche,  F.M.,  and  T.K.  Liu,  "An  Electric  .Model  for  a Ca.ble 
Clamp  on  a Single  Wire  Transmission  Line,"  Report  on 
Contract  F-29601-76-012 5 , AFWL-TR-76-325 , Air  Force  Weaoons 
Laboratory,  .Kirtland  AFB,  NM,  July  1976  . 

2.  Lam,  John,  "Equivalent  Lumped  Parameters  for  a Bend  in  a 

Two  Wire  Transmission  Line : Part  I.  Inductance;  Part  II. 

Capacitance,"  Report  on  Contract  F-29601-76-0125 , AFWL-TR- 
77-5,  Air  Force  Weapons  Laboratory,  Kirtland  AFB,  N.M, 

December  1976. 

3.  Liu,  T.K.,  "Electromagnetic  Coupling  between  Multiconductor 
Transmission  Lines,"  Report  on  Contract  F-29601-76-0125, 
AFWL-TR-7  6- 3 3 3 , .Air  Force  Weapons  Laboratory,  Kirtland 
AFB,  N.M,  December  1976  . 

4.  Tesche,  F..M.,  M.A.  Morgan,  and  B.A.  Fishbine,  "Evaluation 

of  Present  Internal  EMP  Interaction  Technology:  Description 

of  Needed  Improvements,"  AFWL  EMP  Interaction  Note  264, 

Air  Force  Weapons  Laboratory , Kirtland  AFB,  NM,  October  1975. 
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lumped  inductances  and  capacitances  placed  appronr lately 


along  an  otheinvise  uniform  transmission  line.  This  concept 
is  discussed  in  more  detail  in  ref.  (1). 

One  particular  geometry  that  is  ofren  observed  in 
the  internal  configurations  of  aircraft  cables  is  shown 
in  Figure  1,  where  a single  wire  transmission  line  of 
radius  a and  height  b above  the  ground  plane  passes 
over  a thin  septum  of  height  h.  The  wire  and  the  septum 
are  mutually  perpendicular  and  not  touching.  Such  a 
problem  will  model  a cable  passing  over  a rib  in  an  aircraft 
fuselage  or  wing  root.  By  considering  many  periodically 
spaced  septums,  pass  and  stop  band  characteristics  can 
be  determined  for  the  line,  as  outlined  in  ref.  (2). 

In  the  treatment  of  this  problem,  it  will  be 
assumed  that  the  septum  thickness  is  very  small.  This 
implies  that  the  major  effect  on  the  transmission  line 
behavior  will  be  due  to  a capacitive  term  in  the  lumped 
parameter  representation  of  the  obstacle.  The  equivalent 
circuit  of  the  septum,  discontinuity  can  then  be  repre- 
sented, as  shown  in  Figure  2. 

This  paper  describes  in  detail  the  calculation  of 
this  equivalent  capacitance  of  the  septum  and  presents 
the  results  of  a parametric  study  of  this  canonical 
problem. 
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SECTION  II 


FORMULATION  OF  COUPLED  INTEGRAL  EQUATIONS 

The  equivalent  capacitance  of  the  septum  discontinuity 
shown  in  Figure  2 may  be  obtained  by  solving  Laplace's 
equation,  subject  to  certain  boundary  conditions.  If  a 
Green's  function  approach  is  used,  the  problem  may  be 
reduced  to  the  solution  of  a set  of  coupled  integral  equations 
for  the  unknown  excess  charge  distribution  on  the  wire 
and  the  c.harge  distribtuion  on  the  septum. 

The  three-dimensional  free-space  Green's  function 
is  given  by 

G^ir/r')  = 

where  is  the  permittivity  of  free-space,  r is  the 

radius  vector  to  a potential  point,  and  r'  the  radius 
vector  to  a charge  point.  Using  image  theory,  the  perfectly 
conducting  ground  plane  may  be  incorporated  in  the  half- 
space Green's  function  G,  given  by 

G(r/r')  =G^ir/r')  - G^(r/rV  ) (2) 

where  r/  is  the  radius  vector  to  an  image  charge  point. 

The  superposition  principle  now  shows  that 


t 


J"  G(r^/rp  + f 02(?^)  G(rj^/r^)  dS^  = V, 


J G(r2/rj^)  dS|  + J~  02(^2'^  G{r^/rp  dS^  = 0,  ?2^‘”’2 


where  is  the  surface  charge  density  on  the  wire, 

a -y  is  the  charge  distribution  on  the  septuir, , V is  the  poten- 
tial on  the  wire  (with  respect  to  the  ground  plane) , 7^  and 
72  correspond  to  the  surface  of  the  wire  and  diaphragm 
respectively.  The  coupled  integral  Equations  (3)  state 
that  the  potential  on  the  wire  or  the  diaphragm  is  pro- 
duced by  the  charge  distribtuions  ^nd  32- 

The  charge  distribution  contains  two  parts  : 

a uniform  charge  distribution  a i^  the  absence  of  the 

o 

septum  and  an  excess  charge  distribution  ca^  caused  by 
the  septum  discontinuity,  i.e.. 


j , = c + 6a  , 
1 o 1 


This  allows  Equation  (3)  to  be  written  as 


»l 


r-7»^ - 


J dS^  ^ J G 


12  dS'  V 


■ 1 


• 2 


^if-i 


(5a) 


22 


dS. 


(r|)G2idS^ 


f 5b) 


where  G^^^,  = G(r^/rj!^),  Note  that  the  integral  on  the  right  side 
of  Equation  (5a)  defines  the  potential  V on  the  uniform  line, 
and  the  right-hand  side  of  Equation  (5b)  is  the  potential 
i(r2/r^)  at  r2  due  to  the  uniform  line  at  r|  . Equations 
(5a)  and  (5b)  can  be  reexpressed  as: 


J 5ai(r-)  Gii  dS^  +J  dS’  = 0,  ^i^^i 

' 1 ~2 

Jj>7^(rp  G21  dS|  +f  7,  (ry,  G22  dS^  = 

'’1  "2 


(6a) 


(6b) 


where  is  known  exactly  by  suitable  conformal  trans- 

formation (cf.  ref.  5).  Equations  (6a)  and  (6b}  are  the  desired 
coupled  integral  equations  where  s^i  and  ^2  soucht. 

!'Jote  that  these  equations  are  exact. 

5.  Silvester,  P.,  .•^odern  Electromagnetic  Fiells,  '^rentice- 
Hall,  Inc.,  Enqlewood  Cliffs,  'I.J.,  l'^<^8. 
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SECTION  III 


APPROXIMATE  SOLUTION  OF  THE  COUPLED  INTEGRAL  EQUATION 

A rigorous  solution  of  Equations  (6a)  and  (6b)  involves 
integrals  of  the  form 

3=  TO  2~ 


where  6 and  6 ' are  standard  polar  cylindrical  coordinates 
cn  the  wire.  Equation  (7)  is  obtained  by  assuming  that 

OO 

5cJ2^(^;,e')  = ^ cos(n9')  (8) 

n=0 

This  is  equivalent  to  writing  the  Fourier  series  expansion 

for  oj-i  with  coefficients  a that  are  functions  of 

Note  that  the  integral  over  [0,2t]  has  a logarithmic 

singularity  for  and  6=0';  it  cannot  be  integrated 

exactly  for  all  values  of  n.  An  approximate  solution, 

which  is  equivalent  to  truncating  the  Fourier  series 

expansion  at  n=0,  is  used.  This  is  often  referred  to  in  t.he 

literature  as  the  "thin-wire  approximation",  and  is  an 

adequate  approximation  provided  that  the  wire  radius  is  small 

compared  to  other  dimensions  of  the  problem.  A second 

approximation  is  introduced  by  assu.ming  that  t.he  charge 

distribution  on  the  seotum  is  uniform  in  the  y direction. 

* ■'  o 
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Mote  that  the  charge  distribution  on  the  seotum  has  a 
square  root  singularity  at  the  edge  of  the  septum;  the  last 
approximation  may  be  thus  regarded  as  a first  order  approxi- 
mation. The  variational  properties  of  the  capacitance, 
however,  lead  to  a second  order  approximation  in  the  capaci- 
tance results.  .Here,  (x  , y , z„j  denotes  the  actual 

o o o 

dimension  in  cartesian  coordinates. 

Applying  the  thin-wire  approxi.mation  to  Equation  (6) 

we  have 


J io^(rp  dS^  - WT'  / ^ll*^o'^"o^  *^^0 


(9a) 


* 1 


» 


(9b) 


• 1 


where  ot^(x^)  is  a linear  charge  density  related  to 
i5ai(rj^)  via  the  relation 


6t^(x')  = 2ra  3aj_(rp 


(10) 


and  the  kernels  K-]^^  and  K2^  ate  given  by 


(11) 
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V*;’  * ♦ "o  * 


Now  assuming  that 


C2(rp  = f(y<;)  T2(^i) 


where  ^ linear  charge  density  and  ^ (Yq)  is 

uniform  along  the  y'  direction  and  has  dimensions  of 

^ o 

Coul/m;  Equation  (6)  now  shows  that 


■■^2  ^^2^  "^12 


as  n 


- oo  h 


where  the  kernels  K, , and  K_^  are  given  by 


>'l2'V'’i-"i' 


+ (b+v’)^  + (2')^ 

r\  •*  O 


(17) 


Note  that  f(v')  has  been  assumed  constant  and  has  been 
' o 


incorporated  in 


T^.  Tj  aimensicns  of  ooui/m-. 


Equations  (9)  through  (17)  allow  the  coupled  intecral 
equations  (6)  to  be  written  as 


JO 

/ 


00  n 

/ I 


K,  - (X  /•/'  . z'  )dy' 
12  o'^  -o  o o 


= 0, 


< x^  <00  (IE a) 


00 


> ‘'21<'’o' V’'i>  * 


.00  n 

/ T*  (Z')dz'  r K-t  (v^,  z /y  ' , Z' ) dv  ' 
J 2 o o J 22  *o  c *o  o -Q 


= v(y^,  Zq) 


0 1 : h,  — < < » 


(13b) 


For  computational  purposes  it  is  convenient  to 
introduce  dimensional  variables  and  functions  given  by 
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= b^X,y,Z/X-  ,y  ,2-) 


■‘(z;)  = Vb*-  r,(zM 


(19) 


" f'/o' 


L.  .(r/r-),  i,j  = 1,2 


With  this  choice  of  variables  and  functions,  the 
coupled  integral  equations  (18)  take  the  forirt 

® “ h/b 

J :(x')  L^j^(xZx')  dx'  + J n ( z'  ) dz’  J x/y  ’ , z' ) dy  ’ = 0 


— oc  < X < 00 


(20a) 


» 


^ h/b 

(X-)  L23^(y,  2/x’)dx’  + J"  n(z')dz'  J L^^(y,  z/y\z’)dy' 

-»  o 


= -4-r  i (y , z) 


0 f y 1 h/b , - » < z 


(20b) 


Note  that  the  integration  of  L^2  1'22 

(0,h/b)  may  be  performed  exactly  (cf.  ref.  6);  Equation  (20) 
is  thus  reduced  to 

6.  Gradshteyn,  1.3.  and  Ryzhik,  I.M.,  Table  of_ Integrals 
Series  and  Products,  Academic  Press,  New  York,  1365. 
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0 


ni^nL^^{x/z^)dz*  = 


< X < « 


J * J n(3')L*2(y,2/z’)dz' 


= -4'  ;(y,2) 


0 1 y £ h/b , - o < z <? » 


where  L, ^ and  are  given  by 


L*  (X/Z-)  =;.n  •*•  (1-h/b)^  ^ (z’)^  - (I  - h/b) 

12  ^ /-^r - 


>/x^+l+  (z’)“  -1 


_in  Vx^  ->•  (1+h/b)^  + (2')"  -t-  (1  h/b) 


•\/x^  + (2  ' ) ^ + 1 + 


L22(y»z/z')  = 


, „ V (h/b-v) ^ + (z-z')2  + (h/b  - V) 


Vy^  + (z-2’)^  - y 


- 2 n 


(h/b+y)^  + (z-2')^  (h/b  + y) 
fy^  ^ (z-z'K  +y 


-V  «-•***  *»-^'  ^ ^^5 


The  coupled  integral  equation  (21)  may  be  solved 
by  using  the  following  Galerlcin's  approach.  First  inte- 


gratina  Eauation  (21b)  with  respect  to  y aives 


*»  OC 

J -:(x ’)L^^(x/x')dx’  + J 


(x/x')dx’  + / o(z')L^^(x/z’)dz'  = 0 


-•»  < X < ^ 


OC 

I 


o 

* /*  * * 

?(x’)L^,  (z/x')dx'  + / n(z')L,.,  (2/z')dz'  = n(z) 


(24) 


22 


< Z < TO 


If  -k  it 

where  the  )cernels  L2^  and  are  obtained  by  integrating 


22 


and  L22  with  respect  to  y over  (0,h/b).  Similarly,  n is 
obtained  by  integrating  -4:rg  over  (0,h/b).  These  are 
given  by 


L2  ( z/x  ' ) 


L^2  1 


L22  ( z/z ' ) 


h „„  'y/(z-z')^  + (h/b)^  + (h/b) 

b 

(h/b) 


V(z-z')  ^ + (h/b)  ^ - 


+ 2h  V(z-z')  " (h/b)  -*•  (h/b) 

b 


V(z-z’)^  + (2h/b)^  + (2h/b) 


* 3 !z-z'  I + V(z-z'  ) ^ + (2h/b)"^  - 4 yjiz-z')^  + (h/b)  “ 

(26) 
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•a"*  , 


"i  (z ) 


In  (a/2b) 


(1+2^) " 

[ (1-h/b) ^ + z^][ (1+h/b) ^ + z^] 


, h (1-h/b) ^ 

+ _ tn  j 

° (1+h/b)  ^ + z*" 


2tan  ^ (1/z) 


- 2.  [tan.-^  m^)- 


where,  again,  use  was  made  of  table  of  integrals  (ref.  6). 

★ ★ 

Note  that  L^2  (z/z')  has  a logarithmic  singularity 
-In  'z-z'  and  (x/x')  has  a sharp  peak  at  x=x'. 

For  computational  purposes,  these  may  be  treated  as  follows 
First  write 


L22  (z/z')  = 


- i z-z  ' ' + F (z/z  ' ) 

b 


where  F(z/z')  is  continuous  throughout  the  interval 

Next  integrate  by  parts  the  integral  in  ecuation  (24) 
which  contains  L, . (x/x')  . Equation  (24)  then  becomes 


*3 

f (x' 

J 


- (x/x')  dx'  + J n(z')  L^2  (-”/z')  dz'  = 0 

-JO  < X ‘ » 


(x')  L21  (z/x')  dx'  + I n(z')  F(z/z')  dz ' 


JO 


n 2-2 ' n { 2 ’ ) d2  * = -{z 


—X  < 2 ^ ^ 


- 


where 


* , , W (x'  -X)  ^ + (a/b)  ^ (x'-x) 

(x/x')  = >n  -■  ■ 

"V  (x ' -x) ^ + 4 + (x'-x) 

Equation  (29)  is  noted  to  have  a smoother  inteqrand  than 
Equation  (24),  wj'.irh  is  more  amenabl®  to  numerical  computations 
the  integrand  containing  the  loaarithm  will  be  evaluated  in 
closed  form  in  Section  V.  The  second  stage  in  the  use  of 
Galerlcin's  approach  for  obtaining  the  solution  to  Equation  (29) 
requires  knowledge  of  the  properties  of  Hermite  polynomials 
and  the  Gauss-Hermi te  formula.  These  will  be  discussed  in 
the  following  section. 
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SECTION  IV 


HER.MITE  POLYNOMIALS 

The  Hermite  polynomials  Hj(C)  are  defined  over  the 
infinite  interval  and  satisfy  the  orthogonality 

conditions 


CO 


(U  H . (O  dC  = /it  2"-  i!  0 . . 
1 i: 


where  6^j  is  the  Kronecker  delta. 

The  first  few  Hermite  polynomials  are; 


H^(5)  = 1,  H-  (C)  = 2Z.  H^(t)  = A-/  - 2 


and  the  recurrence  relation  is 


Hn^l(C)  - 2?  + 2n  H^_3_(C)  = 0 


These  arise  in  integration  over  and  the  Gauss- 

Her.mite  formula  for  approximating  the  integral  is  given  by 


e"^  f(C)  dC 


= E '"k 


where  is  the  k^  zero  of  and  the  weiahts 

Wj^  is  given  by 


_ 2’^"^  n! 

^ ^ 2 2 


The  weights  w,  and  abscissas  %,  are  given,  for 

K K 

example,  by  Abramowitz  and  Stegun  (ref.  7). 


7.  Abrcimowitz,  M.  , and  J.A.  Stegun,  Handbook  of  Mathematical 
Functions , New  York:  Dover  Publications,  1964. 
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SECTION  V 


NUMERICAL  SOLUTIONS  OF  THE  COUPLED  INTEGRAL  EQUATIONS 
The  coupled  integral  equation  (29)  is  reduced  to  a 
set  of  algebraic  equations  by  using  a Galerkin's  approach 
with  Heraiite  polynomials  as  base  functions.  Assuming  that 


y a H (x') 
/ ■ n n 


n ( z ' ) = e * ^ ' V b H ( z ' ) 
/ . m m 

m=0 


where  a^^,  b^  are  coefficients  to  be  determined,  and  I!,,  (•) 
is  Hermite  polynomial  of  order  Z . Further,  note  that 


■00 


S(z/z')  dz'  (36) 


where  S ( z/ z ' ) = ^ ( z ' - z ) ( ;,n  |Z'-Z'  - 1].  T.his  is  obtained 
by  integration  by  parts  and  the  recurrence  relation  for 
Hermite  polynomials.  The  i.ntegral  equation  (29)  now  takes 
the  form 
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L f 


(x' ) ^ , 


Hn^l(x')  L^^(x/x')dx' 


n = 0 


I / 

m=0 


, ( 2 ' ) 2 * 

b / e ^ H (x')  L,  _{x/z')dz  = 0 ^ x < 

m I m 12 


E / 

n=0 


e ^ H^(x')  L22^  (z/x' )dx' 


M » /» 

Z ^ 0 

m=0 


I ^F(z/z')H^(z 


’ )jdz' 


= r-(  z)  -CO  < 2 < CO 


(37) 


The  integrals  nay  be  evaluated  approximately  by  using 
the  Gauss-Hermi te  quadrature  (33) ; the  result  is 


y a g^^Nx)  ^ y b^  g'^Nx)  = 0 

^ n ’n  ni  ’m 


— CO  < X < =° 


n=0 


m=0 


(38) 


y a (z)  + y b g*'^^  (z)  = f (z) 

^ n ^n  / ■ m 


-CO  < 2 < CO 


n=0 


m=0 


where 
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(39a) 


M, 


"n 


k=0 


L^^(x/x,;)  H^(x^) 


" m 


k=0 


W, 


g'3)(z)  = 


-"E 


^21 


k=0 


M 

1 


= E Wkis'"''"k'  'Vi<'i;>  *'“<"''"k>  v^i> 


k=0 


and  M^  is  the  order  of  this  quadrature. 

Applying  Gauss-Hermite  formula  (33)  once  more  to 
Equations  (38)  results 


E 

n=0 


a . y b T<2) 


n jn 


E 

m=0 


m 


= 0,  j = 0,1,2,,.,; 


N 

E 

n=0 


M 


n Z n 


E 

m=0 


+ > b 

m ?,n 


= d, , z = 0,1,2, .. .N 


23 


(39b) 


(39c) 


(39d) 


(40) 


Equation  (40)  is  a set  of  linear  equations  for  the 


determinations  of  the  unknown  coefficients  a , b , 

n n 

n = 0,1,2, ...N,  which  may  be  solved  on  a digital  computer 

by  means  of  standard  matrix  inversion  routine. 

Once  the  coefficients  a , b^  are  determined,  the 

n m ’ 

capacitance  o':'  the  discontinuity  may  be  obtained  by 

integrating  the  -jxcess  charae  density  on  the  wire. 

The  orthogonality  conditions  of  Hermite  polynomials  may 
be  used  to  perform  this  integration  exactly;  the  result  is 


Notice,  therefore,  that  the  capacitance  depends  only 

on  the  first  term  of  the  Hermite  polynomials  expansion  for 
the  excess  charge  distribution.  In  fact,  it  can  be  shown 
that  t.he  capacitance  of  the  discontinuity , as  given  by 
Equation  (41) , is  stationary  with  respect  to  arbitrary  small 
variations  in  the  functional  form  of  the  excess  charae 
distribution  and  is  lower  bound  (c.f..  Section  IV  of  ref.  '?). 


8.  S.  Coen  and  G.M.L.  Gladwell,  "A  Legendre  approximation 
method  for  circular  microstrip  disk  problem,"  IEEE 
Transactions  on  Microwave  Theory  and  Techniqueil  Vol  .MTT-2S, 
NcT  T',  January  1971.  ^ 
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SECTION  VI 


NUMERICAL  RESULTS 


f 


No  exact  or  approximate  results  to  the  present 
problem  are  known  to  the  authors.  It  is  intended,  there- 
fore, that  the  results  for  the  equivalent  capacitance  of  the 
septum,  obtained  from  the  present  analysis,  be  compared 
with  experimental  results  in  the  future  (ref.  9) . 

Figures  3 through  5 present  dimensionless  excess  charge 
density  on  the  wire,  obtained  from  the  present  analysis, 
each  for  a particular  h/b  and  for  a/b  = 0.00.,  0.01,  0.1. 
Note  that  the  excess  charge  density  is  an  even  function  of 
X , thus  only  the  positive  range  of  is  shown.  These 

results  have  been  obtained  with  N = 10  in  Equation  (35) , 
and  the  order  of  the  Gauss-Hermite  quadrature  is  M^  = 22 
in  equation  (39).  Note  that  the  curves  disolace  kinks  around 
X = 1.6  . This  phenomenon  is  more  evident  for  larcer  values 
of  h/b  and  larger  values  of  a/b.  It  is  attributable  to 
the  fact  that  only  a finite  order  of  the  Gauss-Hermite  quadra- 
ture is  taken.  One  would  obtain  smoother  curves  if  hicher 
order  quadratures  were  used. 

Figure  6 presents  the  equivalent  capacitance  of  the 
septum  C^  , normalized  to  the  capacitance  of  the  uniform  line 

9.  Dr.  Larry  Scott,  Mission  Research  Corporation,  Albuquerque, 
New  Mexico,  private  communication,  December  1976. 
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•xp"' 


(in  the  absence  of  the  septum  discontinuity)  per  unit  length 
times  the  height  of  the  thin  wire  from  the  ground  plane,  as  a 
function  of  h/b  for  a range  of  a/b.  The  same  plot  is 
presented  in  Figure  7 with  h/b  as  a parameter  and  a/b  as 
variables  in  logarithmic  scale. 


Figure  3 


Distribution  of  norinalized  excess  charge 
for  h/b  = 0.1.  Mote  that  is  an 

even  function  of  x. 


27 


0 1 2 3 


X 

Figure  4.  Distribution  of  normalized  excess 
charge  for  h/b  = 0.4. 
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